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Popular Summary 

The summer monsoon is the main source of fresh water supply to the vast population of Asia. 
Recently there have been a growing number of numerical modeling studies indicating significant 
impacts of dust and black carbon aerosols in altering the Asian summer monsoon rainfall 
distribution and large scale circulation. Yet studies of long-term rainfall data over Asia have so 
far shown no clear signals implicating impacts of aerosols. This is partly due to the diverse 
possible forcing of the monsoon climate, from sea surface temperature anomalies, to Eurasian 
snow cover to global warming, which tend to confound the aerosol impact signals. In this paper, 
we used a fingerprinting approach, based on a recent hypothesis of Asian monsoon system 
response to dust and black carbon forcing, to analyze Indian rainfall and related observations to 
detect signals of aerosol impacts. We find that increased absorbing aerosols in the Indo- 
Gangetic Plain in recent decades may have led to long-term warming of the upper troposphere 
over northern India and the Tibetan Plateau, enhanced rainfall in northern India and the 
Himalayas foothill regions in the early part (May-June) of the monsoon season, followed by 
diminished rainfall over central and southern Indian in the latter part (July-August) of the 
monsoon season. These signals, which are consistent with current theories of atmospheric 
heating and solar dimming by aerosol and induced cloudiness in modulating the Indian monsoon, 
would have been masked by conventional method of using all-India rainfall averaged over the 


entire monsoon season. 
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Abstract 


2 

3 In this paper, we present corroborative observational evidences from satellites, in-situ 

4 observations, and re-analysis data showing possible impacts of absorbing aerosols (black carbon 

5 and dust) on subseasonal and regional summer monsoon rain Gill over India. We find that 

6 increased absorbing aerosols in the Indo-Gangetic Plain in recent decades may have led to long- 

7 term warming of the upper troposphere over northern. India and the Tibetan Plateau, enhanced 

8 rainfall in northern India and the Himalayas foothill regions in the early pail (May-June) of the 

9 monsoon season, followed by diminished rainfall over central and southern Indian in the latter 

10 part (My- August) of the monsoon season These signals, which are consistent with current 

1 1 theories of atmospheric heating and solar dimming by aerosol and induced cloudiness in 

12 modulating the Indian monsoon, would have been masked by conventional method of using all- 

1 3 India rainfall averaged over the entiftinonsoon season. 
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2 1. Introduction 

3 The detection and attribution of long-term trends in regional rainfall over the Indian subcontinent 

4 are essential for assessing the consequences of regional climate change, and societal impacts on 

5 approximately 1.2 billion people, about one-sixth of the world population. There is now a 

6 growing body of evidence, mostly from climate model experiments, suggesting that increased 

7 loading of absorbing aerosols (i.e., black carbon, dust, and their mixture) over the India 

8 subcontinent may have strong impacts on monsoon rainfall and circulation [Menon et al., 2002; 

9 Ramanathan et al., 2005; Lau and Kim, 2006; Lau et al., 2008, 2009; Meehl et al., 2008; 

10 Bollasina et al., 2008; Collier and Zhang, 2009, Krishnarnurti et al., 2009; Wang et al., 2009], 

11 While the long-term reduction in surface solar radiation associated with increasing aerosol 

12 loading due to industrial development and population growth in India, has been documented in 

13 many previous studies [Ramanathan et al., 2005; Kumari e/a/, 2007), results on long-term trends 

14 in Indian monsoon rainfall have been varied, depending on the length of the data period, and 

15 whether all-India or regional rainfall data were used [Ramanathan et al., 2005; Joseph and 

16 Simon, 2005: Goswami et al., 2006; IPCC, 2007; Kuniar et al., 2006], More recently, Dash et 

17 al. (2009) has subdivided monsoon rainfall into different regions and seasons, and found 

18 different trends in different regions and different seasons. Most important, all previous studies 

1 9 of Indian rainfall trend are purely statistical in nature with no attempt at physical attributions. In 

20 this study, we use a fingerprinting approach to detect regional and subseasonal bends in Indian 

21 monsoon rainfall, with attribution to possible impacts of absorbing aerosols. 

22 

23 Aerosol may reduce Indian monsoon rainfall by cooling the surface via a reduction of solar 

24 radiation reaching the earth’s surface - the so-called solar- dimming effect (SDM), thereby 
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1 reducing the north-south thermal contrast, resulting in a spin-down of the local monsoon 

2 meridional circulation [Ramanathan et al., 2005]. On the other hand, Lau et al. [2006, 2008] 

3 proposed the “elevated heat pump” (EHP) hypothesis which posits that accumulation of dust and 

4 black carbon over the Indo-Gangetic Plain and the foothills of the Himalayas during the pre- 

5 monsoon season, may induce enhanced warming in the middle and upper troposphere and 

6 subsequently lead to increased rainfall over northern India in late spring and early summer. As 

7 a secondary effect, the increased rainfall and associated changes in cloudiness and water cycle in 

8 the early monsoon may be consequential to reduced monsoon rain over all India in late summer. 

9 They argued that the EHP mechanism may be important to monsoon rainfall not only on 

10 seasonal -to-interannual timescales, but also on long-term (decadal and longer) timescales. 

1 1 Recently, observational studies have showed that a long-term spring warming trend in the upper 

12 troposphere over northwestern India and Pakistan are consistent with the EHP hypothesis 

13 [Gautam et al., 2009; Prasad et al., 2009], However, the passible regional re-distribution of 

14 rainfall over the Indian subcontinent during different phases of the summer monsoon season and 

1 5 relationships to SDM and/or EHP effects are still unknown. 

1 7 2. Data and Approach 

18 The datasets used lor this study include long-term visibility from New- Delhi from the global 

19 surface summary of day data archived at the National Climate Data Center (NCDC), rainfall 

20 from the Global Precipitation Climatology Project (GPCP) for the period 1 979-2007, and from 

21 the Indian regional rainfall data from Indian Institute of Tropical Meteorology (IITM), for the 

22 period 1955-2006. The GPCP rainfall is satellite measurement combined with selected rain 

23 gauge data over land (Huffman et al. 1997). Also used in the study is the Aerosol Index (Al) 
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from Total Ozone Mapping Spectrometer (TOMS), and from the Ozone Measuring Instrument 
(OMI) for absorbing aerosols. The AI is an index that detects the presence of UV-absorbing 
aerosols such as dust and soot. It is defined as the difference between the observations and 
model calculations of absorbing and non-absorbing spectral radiance ratios. . Positive values of 
Aerosol Index generally represent absorbing aerosols (dust and smoke) while small or negative 


values represent non-absorbing aerosols [Hsu et al., 1999], Also used are upper tropospheric 
temperature data from the Microwave Sounding Unit (MSU), specifically temperature profile 
from MSU channel TTS (Temperature Troposphere/Stratosphere) which corresponds to the 
average temperature between 5 to 15 km [Mears and Wentz, 2009] for 1987-2008. Winds from 
the National Center for Environmental Prediction (NCEP) reanalvsis [Kanamitsu et al, 2002] 
will be used to identify associated circulation pattcrnsfiThese data are used collectively to 
establish the fingerprints of aerosol impacts on the Indian monsoon environment. To avoid 
confusion, here we define a trend as a monotonic increase or decrease within the data period as 
represented by the linear regression for the given data period. This definition leaves the 
possibility that a trend unveiled he|& a multi-decadal variation much longer 
than the data record we use. The investigation of such long-term variation is outside the scope of 
this study. 
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3. Results 




3. 1 Increasing absorbing aerosols over the IGP 

Growing population and industrialization in recent decades have led to the rapid increase of 
atmospheric loading of black carbon over the Indo-Gangetic Plain (IGP) - a densely populated 
basin in northern India bounded by the Tibetan Plateau to the north and the Vindhya Range and 
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3. Chota Nagpur Plateau to the south. In addition, during the months of April through June, the 

2 !GP is buffeted by dusty air masses transported by winds from the Thar Desert in northwestern 

3 India and from deserts in Pakistan, Afghanistan, and as far as the Arabian Peninsula, and North 

4 Africa [Fan and Kim, 2006; Lau et al ., 2009], During the monsoon period (Jun-July), dust 

5 particles continued to be transported into the Indo-Gangetic Plain from the Middle East deserts, 

6 across the Arabian Sea. Dust aerosols transported into the IGP may sweep up the finer black 

7 carbon aerosols, which form a coat over the larger dust particles, making them more absorbing 

8 [Eck et al. , 2008; Ramana et al, 2005], The steady increase in atmospheric loading of 

9 absorbing aerosols over the IGP can be seen in the apparent positive trend of Al over the region, 

10 from TOMS and from OMI (Fig. la). Because different instruments are used for different 

1 1 periods, the statistical significance of the Al trend cannot be established in its own right, and can 

12 only be ascertained based on consistency with independent datasets, as well as physically 

13 consistent with model predicted scenario [Lau et al., 2006]. An alternate measure of 

14 atmospheric loading of aerosols is atmospheric visibility, which is directly a function of solar 

15 attenuation by aerosols free of interference from clouds, and is inversely proportional to the 

16 aerosol optical thickness, i.ei', the higher the aerosol loading, the lower the visibility. Long-term 

17 record of visibility is generally available at meteorological stations as a part of the weather 

18 record, A recent study [Wang et aL, 2009a] has established that visibility is a useful proxy for 

1 9 long-term atmospheric aerosol loading both regionally and globally. As a supporting evidence of 

20 the increasing Al trend, the visibility for New Delhi in Fig. la (with axis reversed) also shows a 

21 highly significant trend with greater than 99% confidence level (c.I.) of reducing visibility ( 

22 increasing aerosol loading), that matches well with TOMS Al trend. The long-term reduction in 
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1 visibility is associated with increased attenuation of solar radiation by aerosol, and can be 

2 considered as an indication of the SDM effect [Wang et al., 2009a], 

0 

4 3.2 Large-scale trend patterns 

5 To identify the spatial fingerprints of aerosol forcing on the large-scale environment, we have 

6 computed the linear regressions of the visibility over New Delhi, with a number of large-scale 

7 environmental fields over the Indian subcontinent and adjacent areas for the overlapping period. 

8 Because of the dominant trend in the visibility data, the regression patterns are mostly a 

9 reflection of the long-term trend signals. The upper tropospheric temperature in April-May (Fig. 

10 lb) shows a pronounced warming anomaly over northern India and the Tibetan Plateau, with the 

11 center over the northwestern IGP, representing the large-scale initial set-up of the large scale 

12 environment associated with increasing aerosols in the IGP, similar to that reported in previous 

1 3 studies [Lau and Kim, Gautam e/ al., 2009]. Associafed.with the warming is a large-scale 

14 anticyclone centered over northwestern India at 200 hPa (Fig. lb), with anomalous easterlies 

1 5 over the Indian subebntimcnt and the Bay of Bcngai|> The warming and anticyclonic flow are 

16 dynamically consistent with heating initiated by solar absorption by deep layers of absorbing 

1 7 aerosols, and maintained by induced convective and circulation feedback [Lau et al ., 2006]. The 

18 regressed rainfall pattern shows increased rainfall in May-June over the Bay of Bengal and east 

19 coast of Myanmar (Fig. Icj The increased rainfall extends northwestward along the foothills of 

20 the Himalayas over the IGP. A second positive precipitation anomaly is found over the Arabian 

21 Sea, off the coast of western India, extending northeastward into northwestern India and the Thar 

22 Desert. The rainfall anomalies are associated with increased low-level westerly flow from the 

23 Arabian Sea to the southern Bay of Bengal, with a clear northward flow' towards the foothills of 
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1 the Himalayas over northern India. The upper and lower level flow indicated an anomalously 

2 strong easterly vertical wind shear. Tire large-scale upper tropospheric warming, circulation and 

3 rainfall features are consistent with those of an enhanced Indian monsoon [Webster and yang, 

4 1992; Lau et al 2000], These features are also reminiscent of the observed characteristic 

5 patterns associated with the EHP based on interannual variability [Lau and Kim, 2006], 

6 suggesting that a similar mechanism may be at work on longer timescales. In subsequent months 

7 (July- August), the monsoon rainfall shows a general decline over most of the subcontinent, 

8 particularly over northeastern and southern India, with the upper tropospheric anticyclone 

9 disappearing and replacing by an upper-air cyclonic circulation over the Bay of Bengal, 

10 signaling a weakened Indian monsoon in late summer (see Fig. SI in Supplementary Material 

1 1 and lurther discussion). 

12 

13 3.3 Regional Indian Rainfall fren^L 1% iglfr 

14 For the regional India rainfall trend analysis, we use the IITM regional rainfall data. The 1ITM 

15 data subdivides all-India rainfall into five regions: rtdtthw'est (NW), west central (WC), central 

16 northeast (CME), northeast fNE)> and peninsular (PN), as shown in Fig. 2a. The total rainfall 

17 increases from south (PN-666 ram) to north (CNE=975 mm). There is a strong contrast in 

18 northern India, with the NT. which ts close to the Bay of Bengal, being the wettest (NE-1394 

19 mm), while the northwest region including the Thar Desert, the driest (NW=-483 ram). From 

20 Fig. 2b, the time series of all-Indian-monsoon rainfall (AIMR), defined as the mean of all the 

21 regional rainfall averaged over the extended monsoon season JJAS, show' large interannual 

22 variability' with a slight negative trend of approximately -7.3 mm/decade. However as shown in 

23 Table 1, this trend is statistically insignificant (<90% c.L). 
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To explore the possible signals in regional rainfall associated with the EHP and SDM effects, the 
trends for monthly and bimonthly mean for May through August have been computed for each 
region. Since the trends for May and June are similar, as are those for July and Aug, only the 
bimonthly mean time series are shown in Fig. 3. The numerical values of the trends and 


corresponding statistical significance are summarized in Table 1. For all regions, and time of the 
season, the time series show large interannual variability relative to the trend. Trends are 
discernable in some, but not all regions. During May-June, all regions, except NE, show positive 
trends with NW and CNE which includes the northwest India and the northern central Himalayas 
foothills adjoining Nepal shows the most substantial increase 1 (> 99% c.i.). These regional 
rainfall anomalies are consistent with the positive GPCP rainfall regression pattern over the same 
region (see Fig. lc). In contrast, during July-August negative trends are found in all regions, 
except the NE. Here, the most pronounced signals are found jp PN (>99% c.l.), and in WC 




NAN 


(>90% c.L), which together cover more than half of the area of the subcontinent. 


.41 

Aill 


Overall, the aforementioned large sd qg pattern. and regional rainfall signals arc consistent with 
the EHP and SDM (including induced cloudiness) effects. During the early monsoon season, 

h.'Al'A 

absorbing aerosols heat the atmosphere, and excite atmospheric feedback via the EHP. which 
beside inducing large-scale'rainfall anomalies over South Asia, shifts rainfall to northwestern 
and northern India and the Himalaya foothills (Lau and Kim, 2006). This is consistent with the 
highly significant positive rainfall trend in NW and in CNE in May-June. While rainfall over 
the subcontinent may be heavy during the monsoon season, they occur intermittently, and often 
with prolonged break periods in relative small regions compared to the size of continent. As a 
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1 result, monsoon rain is not necessarily an efficient wash-out agent for atmospheric aerosols, Eck 

2 et al. [2008] showed that total aerosol optical depth in July- August over Kanpur which is located 

3 at the center of the IGP can be as high as 0.5-6, only about 30 % reduction from its peak values 

4 in May-June. The unabated loading of absorbing aerosols and increased cloudiness from the 

5 enhanced early monsoon may act in concert to amplify the SDM effect, causing the land surface 

6 to cool anomalously leading to diminished land evaporation and a reduction of rainfall in the late 

7 monsoon season [Meehl et al 2008; Bollasina et al., 2008; Collier and Zhang, 2009], The 

8 amplified SDM effect may be reflected in the significant negative trend in PN, and WC regions. 

9 The general pattern of reduced rainfall in late summer (July-August) in the m M record is in 

10 overall agreement with the regression patterns of GPCP rainfall in July -August (see Fig. SI in 

1 1 Supplementary Material). 

12 

13 Notice that because of the above subtleties involved in possible aerosol impacts, the regional 

14 rainfall signals will be rriisked if one uses the A1MR for die extended seasonal means May 

15 through September (MJJAS) or June through September (JJAS), as done in many previous 

16 studies [e.g.. Goswami et at.. 2006: Chung and Ramanathan, 2006; Naidu et al., 2009 j. As 

17 shown in Table 1 (first column), both the MJJAS and the JJAS mean for AIMR (Table 1. first 

18 column) show a slight negative, but insignificant trend. This is because the impact of 

1 9 absorbing aerosols on rainfall in the early monsoon season (May-June) is masked by the relative 

20 larger magnitude in the negative rainfall anomalies in the peak monsoon season. Separating the 

21 AIMR into early and late monsoon does produce a significant positive (>95%) trend in May- 

22 June and a marginally significant (>90%) negative trend in July-August. As the preceding 

23 analyses suggest, the positive trend in AIMR is mainly contributed by the northwestern and 
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1 northern regions, while the negative trend by the central and southern regions. Also can be seen 

2 in Table 1 (last two rows), separating by regions but not by early and late season yields no 

3 significant rainfall trend. Hence, timing and regional spatial distribution are key to the detection 

4 of the aerosol impacts on rainfall over India. Previous studies using traditional seasonal (JJA) or 

5 extended seasonal (JJAS) means have indicated that the Indian monsoon has weakened in the last 

6 50 years (Dash el al. 2009, Joseph and Simon 2005). As Table 1 indicates, these seasonal mean 

7 trends are not significant according to our calculations. Rather, the negative rainfall trends 

8 (weakening monsoon) identified by previous studies identified are likely due to rainfall in the 

9 late monsoon season (July- August), and confined to peninsular India and western central regions. 
10 

I 1 4. Concluding remarks 

12 Our results show that the AIMR has a pronounced long-term (45years) positive trend in May- 

13 June, with contribution mostly from the northern part of the subcontinent. A significant negative 

14 trend in AIMR is found in July- August, with the largest contribution from the central and 

15 southern part of the subcontinent. The extended seasonal means (JJAS or MJJAS) used in most 

16 previous studies of the Indian monsoon tend to mask the regional and intra-seasonal trend 

17 signals, which are essential in identifying the fingerprints of aerosol impacts. The positive 

18 rainfall trend in northwestern an.4 northern central India in May-June, is consistent with the 

19 increasing trend in absorbing aerosol in the Indo-Gangetic Plain, and the EPIP large-scale 

20 fingerprints of a warmer upper troposphere, an anomalous anticyclonic circulation over 

21 northwestern Indian and the Tibetan Plateau, and a strong easterly vertical shear across the 

22 Indian subcontinent, all signaling a strengthened Indian monsoon in May-June. The increased 

23 cloudiness from the enhanced rain in the early monsoon season may amplify the SDM effect, and 
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therefore consequential to the reduced monsoon rainfall in July and August. The above signals 
are masked by the use of conventional seasonal mean, and all-India monsoon rainfall. We note 
also that since aerosol impacts on monsoon rainfall are strongly dependent not only on the 
emission sources, but also on transport and deposition by the large-scale circulation, it is likely 
that tire regional Indian rainfall trends attributed to aerosol effects discussed here may also be 
influenced by interactions of anthropogenic forcing and variability stemming from changes in 
large-scale circulation from climate change. 
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1 

2 

3 

4 Figure Captions 

5 Figure 1 a) Time series of Aerosol Index from TOMS, EPTOMS, and OMI averaged over the 

6 Indo-Gangetic Plain [70-90E, 25-35N], in increasing relative unit in the right ordinate and 

7 visibility (solid black line) over New Delhi in decreasing relative units in the left ordinate, 

8 from 1979-2008. Missing data in AI are indicated by symbol- ‘X\ Spatial patterns of linear 

9 regressions of inverse visibility at New Delhi with b) MSU upper tropospheric temperature 

10 (°C) and 200hPA wind, and c) GPCP rainfall (mm day), and 850 hPa winds. Different size 

1 I green dots indicate different levels of statistical significance in the temperate and rainfall 

1 2 fields. Bold arrows indicate wind signal exceeding the 90% significance level. 

13 Figure 2 a) Geographic distribution of the five rainfall regions and annual rainfall total from the 

14 1ITM rainfall data, and b) time series of all-India ra infalf averaged over June- July- August- 

1 5 September. Solid red line indicates linear regression line. 

16 Figure 3 l ime series and linear regression lines (dotted) of each of the five regional rainfall 

1 7 over India for a) May-Jtme (left panels), and for b) July-August (right panels). Units are in 

1 8 mm/day. 
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Table 1 . Linear trend of regional precipitation over India in various seasons for 196 1-2006. Bold 
face types indicate statistically significant trend exceeding 90% level. (*) and (**) indicate 
significant trends at 95% and 99% level, respectively. 


Precipitation Trend (mm/decade) 


Season 

AIMR 

Northwest 

(NW) 

West Central 
(WO 

Central North 
-East(CNE) „„ 

Northeast 

Mm 

Peninsular 

(PN) 

May-June 

9.8* 

11.6** 

8.2 

!9.5** jjj 

p®4.3 

6.6 

July -August 

-9.6 

-7.3 

-13.7 

-4.0 

j|L4 

-18.2** 

MJJAS 

-4.0 

-2.9 

-15,0 

1 3.0 

ll.S 

-13.0 

JJAS 

-7.3 

-5.2 

-15.7 

S9 ' 

1.5 - 

-14.3 





(a) Visibility(New Delhi) & AI(7QE-90E. 
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Supplementary Material 

In addition to the key pre-conditioning and feedback associated with the build-up of absorbing 
aerosols in May-June, the EHP hypothesis also argued that the increased monsoon rainfall and 
circulation induced by absorbing aerosols in late spring and early summer, may be conducive to 
a reduction of the monsoon rainfall in the peak monsoon period (July-August) because of 
increased SDM effect due the increased cloudiness, and changes in the water cycle associated 
with the strengthening of the early monsoon. Figure S 1 b shows that in association with an 
increase aerosol loading (reduction in visibility) over the Indo-Gangetic Plain in April -May, 
there is air overall reduction in rainfall, most pronounced over the northeastern, northern Bay of 
Bengal, and southwestern part of the subcontinent. The large-scale upper level anticyclone in 
April -May (see Fig. lb) has largely disintegrated, and replaceSjly an upper level cyclone over 
the Bay of Bengal, and mostly upper level westerly fk>w over southern India (Fig. SI a), 
consistent with the rainfall reduction. It is important to note that because of the dominance of 
latent heating in the peak monsoon period, the changes in July-August are likely to stem from a 
subsequent dy namical adjustment of the coupled atmosphere-land-ocean system to the initial 
aerosol induced responses, and may be subject to additional remote forcing such as sea surface 
temperature and heat source and sinks outside the domain. This may explain the far less 
organized rainfall pattern and the upper level and lower level wind circulation compared to those 
in May-June (See Fig. I c). 
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Figure SI Patterns of linear regression of OMI AI (April-May) with a) MSU upper troposphere 
temperature and 200 hPa winds for July-August, and b) GPCP rainfall and 850hPa winds for July- 
August. Different size green dots represent different levbis; of statistical significance. Bold arrows 
represent wind anomalies exceeding the 90% significance level. 
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